Propagation of matter wave solitons in periodic and random nonlinear potentials 
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We study the motion of bright matter wave solitons in nonlinear potentials, produced by periodic 
or random spatial variations of the atomic scattering length. We obtain analytical results for the 
soliton motion, the radiation of matter wave, and the radiative soliton decay in such configurations 
of the Bose- Einstein condensate. The stable regimes of propagation are analyzed. The results are in 
remarkable agreement with the numerical simulations of the Gross-Pitaevskii equation with periodic 
or random spatial variations of the mean field interactions. 
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Introduction. Nonlinear excitations in Bose-Einstein 
condensates (BEG) have attracted a lot of attention re- 
cently. In particular matter wave solitons are interesting 
from the fundamental point of view Jl| . The discovery of 
matter waves solitons in BEG 0, 0, 1^1^ has opened the 
possibility to develope new methods for generating and 
controlling solitons. The investigation of the soliton dy- 
namics in inhomogeneous BEG is of interest, in particu- 
lar with inhomogeneities periodic in time or space. Time 
variations can be achieved with the Feshbach resonance 
management technique and it has been studied in 
New type of solitons can be generated by 
this way, and stabilization of higher-dimensional solitons 
in attractive condensate has been shown. Spatial varia- 
tions have been investigated in the form of a periodic or 
random linear potential. In particular, a periodic opti- 
cal lattice can be produced by counter propagating laser 
beams. Such potentials can be used to control the soliton 
parameters [lO|, or to generate gap bright solitons 5] . 
Propagation in a random linear potential was considered 
from the point of view of the Anderson localization in 
BEG and the observation of the crossover between the 
Anderson localization regime and the nonlinear regime 
[1^ . Dark soliton propagation in a linear random po- 
tential was studied in |l3| . In the recent work the 
properties of stationary localized states in the nonlinear 
optical lattice has been investigated. The problems of 
matter wave soliton propagation, when the mean field 
nonlinearity varies periodically or randomly in space re- 
main open. 

In this Rapid Gommunication wc consider the propaga- 
tion of nonlinear matter wavepackets and waves emission 
in the presence of a new type of inhomogeneities, namely 
under nonlinear periodic or random potential, produced 
by periodic or random variations of the atomic scatter- 
ing length in space. The strength and the sign of the 
interatomic interactions, i.e. the value and sign of the 



atomic scattering length a^, can be varied using the FR 
method. Small variations of an external magnetic field 
near the FR can lead to large variations of as ■ For exam- 
ple if we consider the one-dimensional Bose gas close to 
the magnetic wire, then by small variations of the current 
one can induce spatially random magnetic field fluctua- 
tions. This in turn generates random spatial fluctuations 
of the strength of the interatomic interactions 0, . 
Such variations can be achieved also by the optically in- 
duced FR The Gross-Pitaevskii (GP) equation 
describing such configurations has a periodically or ran- 
domly varying in space mean field nonlinear coefficient. 

The dynamics of quasi-lD nonlinear matter waves is 
described by the GP equation 0, 
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Here is the mean field wavefunction, with J jf/jpfia; = 
N, N is the number of atoms, gio — 2hu!±as, where ll!±_ is 
the transverse oscillator frequency, as{x) = asQ + asif{x) 
is the spatially dependent atomic scattering length. The 
spatial dependence will be assumed to be periodic or ran- 
dom. In dimensionless variables where the distance x is 
measured in units of the healing length ^ = h/ ^Jnoguym, 
with no the peak density, and the time t is measured in 
io = ^/(2c), where c = -^/noguj/m, we obtain the equa- 
tion 
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Matter wave soliton motion in a nonlinear periodic po- 
tential. We consider the case when the incident wave is 
the soliton incoming from the left: 
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where 2i>o and AfiQ are the soliton amplitude and ve- 
locity, respectively. The scattering length has spatial 
periodic modulations, so we have V{x) = VoCOs{Kx), 
Vq — 2asi/aso- For small values of Vb, the solution re- 
sembles the unperturbed soliton with modulated param- 
eters in the early step of the propagation. This is true 
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as long as the radiative emission of matter wave is neg- 
ligible. Using the collective coordinate ansatz, we get 
that the soliton mass is preserved, while the soliton cen- 
ter ( obeys the effective particle equation Qt = —dV/d( 
starting from ^(0) = Co) Ct(0) = 4^o, where the effective 
potential is 



V{C) = -A„i cos{KC), 
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Here Ani is the barrier for the soliton moving in the non- 
linear periodic potential. Note that in comparison with 
the hnear periodic potential V{x)u{x,t) in Eq. |2Il, the 
influence of the nonlinear periodic potential on the soli- 
ton is enhanced. For the broad soliton case K/v ^ 1, the 
enhancement factor a = Ani/Ai is a > 1. For the 

narrow soliton case K/v ^ 1, the enhancement factor is 
a = 4i/2/3. The soliton is moving as a classical particle, 
and it can be trapped at C = 2TTn/K,n — 0,1,2.... In 
the trapped regime the soliton performs the oscillatory 
motion with the small oscillation frequency f2 = \/ AniK. 
The critical velocity for depinning of the soliton starting 
from the minimum of the potential is Vdp = ^/2A^. 

When the soliton width is much larger than the period 
of the nonlinear potential, i.e. K/v 3> 1, the radiation 
emission phenomenon can be divided into two time steps. 
First, the soliton emits a small but quick burst of radia- 
tion which is trapped in the form of soliton shape modu- 
lation. Second, the soliton continues to radiate slowly on 
long time scales ~ Vj^"^, as we shall see below. The first 
step can be described in the " renormalized particle limit" 
|lU]. The dressed solution of Eq. (0) can be searched in 
the form u = Us{x,t)(l + x{x,t)). For v < K, where 
V = 4/io is the soliton velocity, the solution for x is 
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At V K the approximation used for the derivation is 
violated, and the radiation emission should be studied 
more carefully. 

Emission of waves by soliton in a nonlinear periodic 
potential. The soliton propagating under action of pe- 
riodic nonlinear potential emits matter wave radiation. 
When the modulations are weak we can use the pertur- 
bation theory based on the Inverse Scattering Transform 
to calculate the radiation emission |l9l |. By conservation 
of the total mass N and energy H, 
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the soliton parameters during the time interval AT 
are modified according to Aiy — — /i, AT), Afi — 
-G{iy,fi,AT), where 
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Here n(A, AT) is the emitted mass density during the 
time interval AT, and fc = 2A (resp. co = 4A^) is the 
wavenumber (resp. the frequency) of the emitted radia- 
tion. The mass (number of atoms) emitted by the soliton 
is Nrad — J n{X)dX. When Vb <C 1 and the propagation 
times are of order VJj"^ we can calculate the emitted mass 
density and the evolution of the soliton parameters. Dif- 
ferent regimes are possible. 

Regime 1. If the modulation wavenumber K is smaller 
that Vq/^o, then the radiative emission is negligible for 
times of order Vg"^. The soliton parameters are almost 
constant in this regime. 

Regime 2. If the modulation wavenumber K is larger 
than Vq/hq, then the soliton emits a significant amount 
of radiation. The soliton amplitude and velocity satisfy 
the system of ordinary differential equations 



vt = -F{v, /i), fit = -G{u, fi), 
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starting from I'iO) — iyo,fi{0) — /iQ. The functions F and 
G are given by 
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The maximal radiative decay is obtained for K close to 
I \i. By integrating Eq. I^J, we can put into evidence 
that there are two subcases: a). The soliton mass 4j/ 
decays, while the velocity increases or decays slowly, so 
that the soliton stays in the regime K[i > . Eq. can 
be used to describe the long time behavior of the soliton, 
whose mass decays to zero (Figure h). The soliton 
mass decays, but the velocity 4^ decays faster, so 
that the soliton parameters reach the condition K[i = 
at the critical time. This state is a stable equilibrium, 
we recover the regime 1. This shows that we can have a 
stable soliton even in the case K > hq/i'q, at the expense 
of the emission of a small amount of radiation to allow 
the soliton to reach a stable state fFigures 1213(1 . 

Regime 3. If we start from critical initial conditions 
such that hqK = Vq, then there are three sub-cases: a), 
li K < AfiQ, then the soliton is attracted by the regime 
2 and its mass decays to 0. b). If i^T > 4/io, then the 
soliton is attracted by the regime 1 and its parameters 
remain constant, c). If K = 4/io (and thus vq = 2/io), 
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FIG. 1: Soliton amplitude. The initial soliton parameters are 
fc'o = Ato = 1- A periodic modulation 0.1 cos(_R'a;) is applied 
to the nonlinear parameter. We compare the results from full 
numerical simulations of the perturbed GP equation (thin 
dashed lines) with the theoretical predictions of Q (thick 
solid lines). From top to bottom: K = 1.5 (regime 1), K = 
0.5 (regime 3a), K = 1 (regime 2a). The strongest decay is 
achieved for the critical case K — 1 = i/q / fio, <i-s predicted by 
the theory. 
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FIG. 2: Soliton amplitude. The initial soliton parameters 
are vq — 1, fio — 0.5. Left picture, from top to bottom: 
K = 1.5 (regime 1), K = 2.5 (regime 2b). Right picture: 
K = 2 (regime 3c). As predicted by the theory, the critical 
case experiences oscillations but is stable. 



the soliton experiences strong oscillations but its mass 
does not decay (Figure |2Jl. 

Emission of waves by soliton in a nonlinear random 
potential. Let us assume that the function V is the real- 
ization of a random zero-mean stationary process. The 
correlation function is B{x,lc) = {V{x)V{0)), where Ic 
is the correlation length. For times of order B{0,lc)~^ 
the soliton parameters satisfy the deterministic system 
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FIG. 3: Soliton amplitude. The initial soliton parameters are 
vo = 1, /^o = 0.25. Left picture, from top to bottom: K = 2 
(regime 1), K = 8 (regime 2b). Right picture: K = 4 (regime 
3b). As predicted by the theory, the critical case is stable. 



@ where the functions F, G are defined by 
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Here the function c is given by 

((A + M)^ + ^^)((A-A*)^ + ^^) 
cosh(7r(i^2+A2-^2)/(4^j,)) ' 

and the coefficients d and k by: 



(A ~ m)' + 



d{k) — '2 B{x) cos{kx)dx. 
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Note that fc(A) > i^'^/n for all A. Thus, the interaction 
between the soliton and the nonlinear random potential 
depends only the tail of the power spectral density d{k) 
of y for fc > There are two regimes of propagation: 

Regime 1. If //q S> J^Oi then the emitted radiation 
density is concentrated around the wavenumbers ±2/io. 
Besides the system Q can be simplified, and we obtain 
that the velocity of the soliton is almost constant, while 
the mass decays as a power law: 
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In this regime the radiative decay prevents from trans- 
mitting nonlinear wavepackets. The decay time is in- 
versely proportional to the correlation function of the 
disorder and the forth power of the soliton amplitude. 
This means that this type of disorder intensively destroys 
heavy solitons. 

Regime 2. If /io ^ J^Oj then the soliton emits a very 
small amount of broadband radiation, its mass is almost 
constant, while the velocity decays very slowly, typically 
as a logarithm |2(]| . The amount of emitted radiation is 
proportional to diy^ / ^i). In this regime the soliton can 
be transmitted. 

The analysis of the system Q shows that these two 
regimes are attractive, in the sense that, after a tran- 
sient regime, one observes the regime 1 (resp. 2) if iiq/vq 
is above (resp. below) a critical value. We have checked 
these predictions by numerical simulations of the ran- 
domly perturbed GP equation We consider the case 
of a stepwise constant process V . The constant step is 
equal to Ic, and the process V takes random indepen- 
dent values over each elementary interval that are uni- 
formly distributed in [— cr, a]. The power spectral density 
is d{k) = 2(t2[1 — cos(fc?c)]/[fc^Zc]- We compare in Figure 
21 the numerical results and the theoretical predictions 
in a case close to the regime 1 (top figure) and close to 
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FIG. 4: Soliton amplitude. A random modulation with am- 
plitude a and correlation length U is applied to the nonlinear 
parameter. The initial soliton parameters and potential am- 
plitude are = /xo = 1, cr = 0.1 (top) and i/o — 1, fio = 0.5 
a = 0.2 (bottom). We compare the results from full numerical 
simulations of the perturbed GP equation (thin dashed lines) 
with the theoretical predictions of Q (thick solid lines). For 
each figure, we have from top to bottom: Ic — S, h ~ 0.125, 
L = 1. 



the regime 2 (bottom figure). The power law radiative 
decay is noticeable in the top figure, with a decay rate 
that is maximal for Ic — 0.58. The full transmission 
regime (up to a transient regime where the soliton emits 
radiation) can be observed in the bottom figure. 

To estimate the typical values of the parameters, let us 
consider the case of ^Li condensate with the FR a,t B = 
720G. The typical values of the scattering length used 
in the experiments are: at B = 352G is as — — 0.23nm 
and as(300G) « — 0.2nm Thus by changing period- 
ically in space the magnetic field between these values 
we can obtain the nonlinear periodic potential with the 
dimensionless amplitude of modulations Vo ~ 0.26. In 
the trap with uj± — 2% ■ lO'^Hz and uq — 10® cm^^, we 
have ^ « 2/zm, c 2 mm/s, so to ~ 1 ms. A soliton with 
velocity f ~ c travels through the region with modulated 
scattering length with L 0.25mm in the dimension- 
less time t ~ 125. The soliton width is of the order of 
1.5/xm, and if the grating period varies in the interval 
(1, 15)/xm, then both limits Ki/ ^ 1 and Kiy <C 1 are 
covered. The optically induced FR method gives a grat- 
ing period l/im. The critical case K = 1 — i'o/mo 
(Fig. 1) corresponds for these parameters to the grat- 
ing period ~ 6/im and the soliton velocity ~ 8mm/s. 
The random modulations with the deviation strength is 
tts/asi = 0.1 can be achieved by the random distribution 
of the current in wire along the atom chip[l5j. For the 
soliton velocity = 2c, and the initial number of atoms 



in soliton is ~ 10^, we find then the soliton decay time 
Tc ^ 50, which is equal to ~ 0.1s in physical units. 

In conclusion we have investigated the transmission 
of matter wave bright solitons through periodic or ran- 
dom nonlinear potential, generated by periodic or ran- 
dom spatial variations of the atomic scattering length. 
The condition for the emission of matter waves and ra- 
diative soliton decay are obtained. We show that critical 
cases support the stable propagation of bright solitons. 

Acknovifledgements. The authors thank B.B. Baiza- 
kov and M. Salerno for interesting discussions. Partial fi- 
nancial support of the work of F.Kh.A. from the Physics 
Department of Salerno University is acknowledged. 



[1] 
[2] 

[3] 
[4] 

[6: 

[7 

[8 

[9: 

[10 

[11 

[12 
[13 
[14 
[15 
[16 

[17 
[18 



[19 
[20 



V. A. Brazhnyi and V. V. Konotop, Mod. Phys. Lett. B 
18, 627 (2004). 

S. Burger, K. Bongs, S. Dettmer, W. Ertmer, K. Seng- 
stock, A. Sanpera, G. V.. Shlyapnikov, and M. Lewen- 
stein, Phys. Rev. Lett. 83, 5198 (1999). 
K. S. Strecker, G. B. Partridge, A. G. Truscott, and R. 
G Hulet, Nature, 417, 150 (2002). 

L. Khaykovich, F. Schreck, G. Ferrari, T. Bourdel, J. Cu- 
bizoUes, L. D. Carr, Y. Castin, and C. Salomon, Science 
296, 1290 (2002). 

B. Eiermann, Th. Anker, M. Albiez, M. Taglieber, P. 
Treutlein, K.-P. Marzlin, and M. K. Oberthaler, Phys. 
Rev. Lett. 92, 230401 (2004). 

F. Kh. AbduUaev, E. N. Tsoy, B. A. Malomed, and R. 

A. Kraenkel, Phys. Rev. A 68, 053606 (2003). 

F. Kh. AbduUaev and M. Salerno, J. Phys. B 36, 2851 
(2003); F. Kh. AbduUaev, A. Gammal, and L. Tomio, J. 
Phys. B 37, 635 (2004). 

P. G. Kevrekidis, G. Theocharis, D. J. Frantzeskakis, and 

B. A. Malomed, Phys. Rev. Lett. 90, 230401 (2003). 

G. Theocharis, P. Schmelcher, P. G. Kevrekidis, and 
D. J. Frantzeskakis, Phys. Rev. A 72, 033614 (2005); 
arXiv:cond-mat/0509471 

R. Scharf and A. R. Bishop, Phys. Rev. E 47, 1375 
(1993). 

B.B. Baizakov, B.A. Malomed, and M. Salerno, Euro- 
phys.Lett. 63, 642 (2003). 

U. Gavish and Y. Castin, Phys. Rev. Lett. 95, 020401 
(2005). 

N. Bilas and N. Pavloff, Phys. Rev. Lett. 95, 130403 
(2005). 

H. Sakaguchi and B. A. Malomed, Phys. Rev. E 72, 
046610 (2005). 

H. Gimperlein, S. Wessel, J. Schmiedmayer, and L. San- 
tos, Phys. Rev. Lett. 95, 170401 (2005). 
S. Wildermuth, S. Hofferberth, I. Lesanovsky, E. Haller, 
M. Andersson, S. Groth, L Bar- Joseph, R. Folman, and 
J. Schmiedmayer, Nature 435, 440 (2005). 
P. O.Fedichev, Yu. Kagan, G. V. Shlyapnikov, and J. T. 
M. Walraven, Phys. Rev. Lett. 77, 2913 (1996). 
M. Theis, G. Thalhammer, K. Winkler, M. Hellwig, G. 
Ruff, R. Grimm, and J. H. Denschlag, Phys. Rev. Lett. 
93, 123001 (2004). 

V. L Karpman, Phys. Scr. 20, 462 (1990). 

J. Garnier, SIAM J. Appl. Math 58, 1969 (1998). 



